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We have investigated motional heating of laser-cooled 9 Be + ions held in 

radio-frequency (Paul) traps. We have measured heating rates in a variety of 

traps with different geometries, electrode materials, and characteristic sizes. 

The results show that heating is due to electric-field noise from the trap 

electrodes which exerts a stochastic fluctuating force on the ion. The scaling 

of the heating rate with trap size is much stronger than that expected from 

a spatially uniform noise source on the electrodes (such as Johnson noise 

from external circuits), indicating that a microscopic uncorrelated noise source 

on the electrodes (such as fluctuating patch-potential fields) is a more likely 

candidate for the source of heating. 



I. INTRODUCTION 



Cold trapped ions have been proposed as a physical implementation for quantum com- 
putation (QC) [l|], and experiments on one |2|-|5| and two ||[7| ions have demonstrated proof 
of the principle. Work is currently underway to extend these results. In ion trap QC, 
ion-laser interactions prepare, manipulate and entangle atomic states in ways dependent on 
the quantum motional state of the ions. A limiting factor in the fidelity of an operation 
is uncontrolled heating of the motion during manipulations. Heating leads to decoherence 
of the quantum superposition states involved in the computation |8],|9|], and can ultimately 
limit the number of elementary gate operations which can be strung together. Speculations 



have been made about the mechanisms that lead to heating P,|T0|- |14[| , but measurements 
are scarce since the necessary sensitivity can be achieved only through laser cooling to near 
the ground state of motion. Additionally, systematic studies of the dependence of heating 
rate on various trap properties are difficult, since often this requires the construction and 
operation of an entirely new trap apparatus which may have different values of properties 
not under study. Indeed, the data presented here pose several interpretational difficulties 
for this reason. 

Heating of a single trapped ion (or the center-of-mass motion of a collection of trapped 
ions) occurs when noisy electric fields at the position of the ion couple to its charge, giving 
rise to fluctuating forces. If the spectrum of fluctuations overlaps the trap secular motion 
frequency or its micro-motion sidebands, the fluctuating forces can impart significant energy 
to the secular motion of the ion. Here, we express the heating rate as the average number 
of quanta of energy gained by the secular motion in a given time. There are several candi- 
dates worth considering for sources of the noisy fields which give rise to heating. Some of 
these are ||: Johnson noise from the resistance in the trap electrodes or external circuitry 
(the manifestation of thermal electronic noise or black body radiation consistent with the 
boundary conditions imposed by the trap electrode structure), fluctuating patch-potentials 
(due, for example, to randomly oriented domains at the surface of the electrodes or adsorbed 
materials on the electrodes), ambient electric fields from injected electronic noise, fields gen- 
erated by fluctuating currents such as electron currents from field-emitter points on the trap 
electrodes, and collisions with background atoms. Only the first two mechanisms will be 
considered here since the remaining mechanisms (and others) are unlikely contributers || or 
can be eliminated by comparing the measured heating rates of the center-of-mass and differ- 
ential modes of two ions ||. As will be shown below, the Johnson-noise and patch-potential 
mechanisms give rise to heating rates which scale differently with the distance between the 
ion and the trap electrodes. 
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II. TWO MODELS FOR SOURCES OF HEATING 



A. Preliminaries 



The heating rate caused by a fluctuating uniform field can be derived as in Savard, et 
al. [15 1 and agrees with a classical calculation ||. The Hamiltonian for a particle of charge 



q and mass m trapped in an harmonic well subject to a fluctuating, uniform (non- gradient) 
electric field drive e{t) is 

H(t) = H - qe(t)x, (1) 

where H = p 2 /2m + mu^x 2 is the usual, stationary harmonic oscillator Hamiltonian with 
trap frequency u m . From first order perturbation theory, the rate of transition from the 



ground state of the well (|n = 0)) to the first excited state (\n = 1)) is |15 



OC 



rW = i / dre^(e(t)e(t + r))\(0\qx\l)\ 2 . (2) 
n J-oo 

Evaluating the motional matrix element gives 

q 2 

T ^i = — i S E {u m ) (3) 

mnuj m 

where Se{u) = dTe luJmT (e(t)e(t + r)) is the spectral density of electric-field fluctuations. 

For an ion trapped by a combination of (assumed noiseless) static and inhomogeneous 
rf fields (Paul trap) the heating rate can be generalized to JTT| : 

h =^r m { SE(ulm)+ i SE(aT±ul "^)' (4) 

where n is the rate of change of the average thermal occupation number, uo m is now the 
secular frequency of the mode of motion under consideration and Qt is the trap rf-drive 
frequency. The second term on the rhs of Eq. ^ is due to a cross-coupling between the rf 
and noise fields; it will not be present for the axial motion of a linear trap which is confined 
only by static fields. Even for motion confined by rf pondermotive forces, this second term 
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will be negligible in the absence of spurious resonances in Se(uj) or increasing Se{w) (since 
uj^/VLj, ~ 10~ 4 ) and is neglected in what follows ||. 

We differentiate two sources of the noise that gives rise to heating. The first is thermal 
electronic noise in the imperfectly conducting trap electrodes and elsewhere in the trap 
circuitry. Though this source of noise is ultimately microscopic in origin, for our purposes 
here it can be treated adequately by use of lumped circuit models. Thermal noise has been 
considered in the context of ion-trap heating in several places |]8|]il"Hl4l . The second source of 
noise considered here is due to "microscopic" regions of material (small compared to the size 
of the trap electrodes) with fluctuating, discontinuous potentials established, for example, 
at the interface of different materials or crystalline domains. We call this patch-potential 
noise, and its microscopic origin leads to manifestly different heating behaviour from that 
for the thermal electronic noise case. Static patch-potentials are a well-known phenomenon, 
but little is known about the high-frequency (MHz) fluctuating patches which are required 
to account for our observed heating rates |I6|-|IE|] . 



B. Thermal electronic noise 

Heating rates in the case of thermal electronic noise (Johnson noise) can be obtained 
simply through the use of lumped- circuit models, which are justified by the fact that the 
wavelength of the relevant fields (at typical trap secular or drive frequencies) is signifi- 
cantly larger than the size of the trap electrodes. Such an analysis has been carried out 
elsewhere [p| pT| , p!9| , and only the major results will be quoted here. Resistances in the 



trap electrodes and connecting circuits give rise to an electric field noise spectral density 

Se{oj) = 4kBTR(uj)/d 2 where d is the characteristic distance from the trap electrodes to 

the ion, T is the temperature (near room temperature for all of our experiments), k B is 

Boltzmann's constant, and R(u) is the effective (lumped-circuit) resistance between trap 

electrodes. The heating rate is given by 

A = q 2 k B TR(u m ) 
mhtu m d 2 



A numerical estimate of the heating rate for typical trap parameters gives 0.1/s<nR<l/s 
pjJTTJ] , which is significantly slower than our observed rates. As a final note, the lumped 
circuit approach is convenient, but not necessary. In the Appendix, we present a microscopic 
model that is valid for arbitrary ion-electrode distances and reproduces Eq. |5] for all the traps 
considered here (and for all realistic traps where d ^> 5, where 5 is the skin depth of the 
electrode material at the trap secular frequency). 



C. Fluctuating patch-potential noise 

To derive the heating rate for the case of microscopic patch-potentials we use the following 
approximate model. We assume that the trap electrodes form a spherical conducting shell 
of radius a around the ion. Each of the patches is a disc on the inner surface of the sphere 
with radius r^-Ca and electric potential noise V p (uj). Alternatively, each patch is assumed 
to have power noise spectral density Sy(u). The electric field noise at the ion due to a single 
patch is E p (u) = — ?>V p {uj)r 2 / Aa^ in the direction of the patch. There are N « 4Ca 2 /r 2 
such patches distributed over the sphere with coverage C < 1. Averaging over a random 
distribution of patches on the sphere, we find that the power spectral density of the electric 
field at the ion (along a single direction) is 

q f , AT ( dE p (u) \ 2 3CS v (u)t* 



This gives a heating rate 



^ = 3q 2 CrlS v (u; m ) 
lQmhumd 4 ' 



in which the association d ~ a is made. Note the difference in scaling with electrode size 
between Eqs. |5| and 0. The thermal electronic noise model gives a scaling % oc d~ 2 , while 
the patch-potential model gives np oc d~ 4 . In fact, a d~ A dependence also arises from a 
random distribution of fluctuating charges or dipoles. 
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III. MEASUREMENTS 



A. Measuring the heating rate 

To determine the heating rate, we first cool the ion to near the ground state. In suf- 
ficiently strong traps, this is achieved simply by laser cooling with light red-detuned from 
a fast cycling transition (7 rs u m , where 7 is the radiative linewidth of the upper state) 
propagating in a direction such that its k- vector has a component along the direction of the 
mode of interest. In weaker traps, additional sideband Raman cooling is utilized to cool to 
the ground state [[|. Typical starting values of n, the average number of thermal phonons 
in the mode of interest, are between and 2. 

After cooling and optically pumping the ion to its internal ground state (denoted | J,)), 
we drive Raman transitions between atomic and motional levels Tuning the Raman 

difference frequency Auj to the k th motional blue sideband (bsb) at Auj = u Q + ku m drives 
the transition | |}|n) <-> | \)\n + k) where | |), | [) refer to the internal (spin) states of 
the atom that are separated by ojq. The k th red sideband (rsb) at Auj = ujq — koj m drives 
I i)\ n ) *~* I T)l n ~~ The measurement utilizes asymmetry in the strengths of the red and 
blue motional sidebands to extract n. The strengths of the sidebands are defined as the 
probability of making a transition | J,) | f), which depends on the occupation number 
of the motional levels. The strengths are probed by a Raman pulse of duration t tuned to 
either k th sideband. The probability Pj_ of remaining in | {) after probing is measured and 
the strengths Jjf b = 1 — P^ T sb and I^ sh = 1 — P^bsb are extracted. For thermal motional 
states, the strengths of the red and blue sidebands are related by |||8j 



-^fc Sb _ S ^ n2 ^m,rn-kt (8) 

oo 

^ ] Pm sin f2 m ^fc j?Ti t (9) 



m=k 



n x fe oo 



1 + n 



m=0 



4 bsb (io) 



.1 

where Q m +k,m — Qm,m+k is the Rabi frequency of the k th sideband between levels m and 



m + k, and P m = n m /(l + n) m+1 is the probability of the m th level being occupied for 
a thermal distribution of mean number n. The ratio of the sidebands = 7^ sb /J^ sb is 
independent of drive time t and immediately gives the mean occupation number n, 

,- i = (in 

which is valid even if the Lamb-Dicke criterion is not satisfied. In principle, k should be 
chosen to be the positive integer nearest to n in order to maximize sensitivity. In practice 
we use k =1, 2, or 3 in most cases. Note that Eq. [ll| is valid only for thermal states; this 
is adequate since Doppler cooling leaves the motion in a thermal state P,|20], as does any 
cooling to near the ground state. 

In order to determine the heating rate n, delays with no laser interaction are added 
between the cooling cycle and the probing cycle. An example of a data set at a fixed trap 
secular frequency is shown in Fig. [l|. The error bars are determined as follows: The raw 
data of Raman scans over the sidebands (such as is those shown in the insets of Fig. [I]) are 
fit to Gaussians, from which the depths of the sidebands are extracted, with error on the 
parameter estimate calculated assuming normal distribution of the data. The errors from 
the rsb and bsb strengths are propagated through Equation |IT] for an error on n. The error 
bars shown are one-sigma, and include only statistical factors. These errors are incorporated 
in the linear regression to extract n with appropriate error. Many such data sets are taken 
for various types of traps and at different secular frequencies. 

B. The traps 

The measurements of heating rate in this paper extend over a five year period and utilize 
six different traps. The traps are summarized in Table |. The traps are described in the 
references listed; here only a brief discussion is included. The "ring" traps are approximate 
quadrupole configurations consisting of a flat electrode (125 /im thick) with a hole drilled 
through it (the ring) and an independent "fork" electrode (100 fim thick) that forms endcaps 
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on either side of the electrode, centered with the hole, similar to the trap shown in Fig. |2|. 
In trap 1, the ring and endcap electrodes are at the same average potential; in traps 2 and 
3 a static bias field could be added between the fork and ring to change the distribution of 
binding strengths along the three principle axes of the trap. The size of these traps is stated 
as the hole radius, with the endcaps spacing approximately 70% of the hole diameter. For 
the elliptical ring trap (trap 2) the stated size is the radius along the minor axis and the 
aspect ratio is 3:2; the fork tines are parallel to the major axis of the ellipse. Traps 3a and 3b 
were drilled into a single flat electrode with a single graded fork electrode (see Fig. |2|). The 
rings are circular and the size stated is the radius. This was the trap used for the size-scaling 
measurements. The heating in all of the ring traps was measured in a direction in the plane 
of the ring electrode, parallel to the tines of the fork electrode. Traps 4, 5 and 6 are similar 
linear traps with geometry indicated in Fig. [3]. Trap 6 was made slightly larger than traps 4 
and 5 by increasing the space between the two electrode wafers. Heating was measured along 
the axial direction, which has only a static confining potential. The size quoted in Table | 
for the linear traps is the distance between the ion and the nearest electrode. All traps are 
mounted at the end of a coaxial A/4 resonator for rf voltage buildup pi| . Typical resonator 
quality factors are around 500 and rf voltage at the open end is approximately 500 V with 
a few watts of input power. In all traps except for trap 3a and 3b the resonator is inside the 
vacuum chamber with the trap. In traps 3a and 3b, the resonator is outside the chamber, 
with the high-voltage rf applied to the trap through a standard vacuum feedthrough. 

Since we believe that surface effects are an important factor in heating, we cleaned the 
electrode surfaces before using a trap. When trap electrodes were recycled, they were first 
cleaned with HC1 in order to remove the Be coating deposited by the atomic source. For the 
molybdenum traps an electro-polish in phosphoric acid was then used. For the beryllium 
electrodes electro-polishing in a variety of acids was ineffective, so abrasive polishing was 
used. Finally, the traps were rinsed in distilled water followed by methanol. The gold 
electrodes of the linear traps were cleaned with solvents after being evaporatively deposited 
on their alumina substrates. The time of exposure of clean trap electrodes to the atmosphere 
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before the vacuum chamber was evacuated was typically less than one day. The traps were 
then vacuum baked at ~ 350° C for approximately three days. 

C. Data 

Our longest-term heating measurements were made on trap 1. In Figure |] we plot the 
heating rate as a function of date of data acquisition for a fixed trap frequency (11 MHz). 
The heating rate is on the order of 1 quantum per millisecond with a basic trend upwards 
of ~ 1 quantum per millisecond per year. Over this time the electrodes were coated with 
Be from the source ovens, but beyond this, nothing was changed in the vacuum envelope, 
which was closed for this entire period of time. The cause of the increase in heating rate is 
unknown, but may be related to increased Be deposition on the electrodes. Be plating on 
the trap electrodes could be a source of patch-potential noise. 

Figure [| shows heating rates in the linear traps (trap 4, 5 and 6) and the elliptical ring 
trap (trap 2) as a function of trap secular frequency. The frequency dependence of the 
heating rate is expected to scale as S^(u; m ) /cj m (Eq. For example, a trap electrode with 
a flat noise spectrum (Se{u) = constant) will have a heating rate that scales as a;" 1 . The 
actual spectrum of fluctuations is impossible to know a priori, but in principle the data can 
be used to extract a spectrum over a limited frequency range given the model leading to Eq. 
f|. For the three linear traps, the heating rate data are most consistent with a a;" 2 scaling, 
implying Se u;^ 1 . This does not greatly assist in identifying a physical mechanism for the 
heating. For example, pure Johnson noise will have a flat spectrum, low-pass- filtered Johnson 
noise will have a spectrum that decreases with increasing frequency, and the spectrum of 
fluctuations in the patch-potential case is entirely unknown. In addition to the theoretical 
ambiguity, there is evidence in other data sets of different frequency scalings (though they 
are always power-law scalings). This measurement certainly cannot be used to pinpoint a 
heating mechanism; it is presented here only for completeness. 

The data of Figure |5|b provide a first indication of the scaling of the heating rate with 
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trap size. Trap 6 is about 1.3 times larger than trap 5, while its heating rate (at 10 MHz) 
was a factor of 3 slower. This indicates that the dependence of heating rate on trap size 
is stronger than d~ 2 , but is consistent with d~ A . Of course, this comparison is to be taken 
with some caution, since these are two separate traps measured several weeks apart, and 
therefore likely had different microscopic electrode environments. However, a comparison 
is warranted since the traps were identical apart from their sizes. In particular, all the 
associated electronics was the same and the rf drive voltage was very nearly the same. In 
fact, the rf voltage was slightly larger for the measurements on trap 6, which showed the 
lower heating rate. This is important to note because we observe a slight dependence of 
the heating rate on the applied rf trapping voltage. Though we have only a qualitative 
sense of this dependence at present, it seems that heating rates increase with rf voltage, 
up to a point, at which the effect levels off. This rf- voltage dependence is observed along 
directions where the ion is confined both by static fields and by pondermotive fields. It may 
not be unreasonable that the increased rf voltage increases the intensity of the noise source 
(possibly due to an increase of temperature of the electrodes), even when it does not affect 
the trap secular frequency, as in the axial direction of the linear traps. 

Trap 3 was designed to give a controlled measure of the heating rate as a function of trap 
size, while all other parameters were held fixed. The trap electrodes were made from the 
same substrates, the electrodes were subjected to the same pre-use cleaning, the traps were 
in the same vacuum envelope, driven by the same rf electronics (simultaneously) and data 
for both traps were acquired with minimal delay. For direct comparison at the same secular 
frequency in both traps, it was necessary to change the applied rf voltage since ou m oc 1/d 2 . 
(A static bias between ring and endcap can be added, as discussed above, but this was not 
sufficient to measure heating at identical secular frequencies for the same rf drive.) There 
are two data sets to be discussed for this trap, shown in Figure ^. 

In the first set, shown in Figure |]a, we have data points at two different secular frequen- 
cies for the "small trap" (trap 3a) and one point for the "big trap" (trap 3b). The heating 
rates of the small trap are comparable to the heating rates for other traps and show a cu" 1 
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scaling of the heating rate. The single point on the big trap is at a lower secular frequency, 
yet has a much slower heating rate. In fact, if we extrapolate the data from the small trap 
to the same secular frequency (using a;" 1 ), the heating rate is over an order of magnitude 
lower in the big trap. The ratio of the heating rate in the small trap to that of the big trap 
is 20 ± 6. This is a much stronger scaling than that predicted by a Johnson noise heating 
mechanism (Eq. [| predicts a d 2 ~ 4.8 scaling), but is consistent with the scaling in the 
patch-potential case (Eq. |7| predicts a d A ~ 23 scaling). When these data are used to predict 
an exponent for the size-scaling, the result is <i 3,8±0 ' 6 . 

For the second data set, shown in Figure ||b, the trap was removed from the vacuum 
enclosure, given the usual cleaning (as discussed above), and replaced for the measurements. 
In this data set, the trap behaved quite differently from all other traps, with heating rates 
significantly below those of other traps. Also, Se must have been a strong function of u> for 
this trap since the scaling with trap frequency was rather pronounced. The scaling with size 
was also strong: the heating rate was 16,000 times smaller in the big trap. When these data 
are used to predict an exponent for the size-scaling, the result is d 12±2 . Needless to say, it is 
difficult to draw general conclusions from the data for this particular trap, but the difference 
in heating rates between the two traps seems to strongly indicate, again, that Johnson noise 
is not the source of the heating. We cannot be sure why this trap had such anomalous 
heating behaviour, but we speculate that it is due to a less-than-usual deposition of Be 
on the trap electrodes prior to the measurements, because the trap loaded with minimal 
exposure to the Be source atomic beam. 

At this point it is useful to compare the present results to heating rates in other ex- 
periments. There are two other measurements. The first was done with 198 Hg + [f[0[| . For 
that experiment u; m /27r « 3 MHz and d ~ 450 /im and the heating rate was 0.006/ms. 
Accounting for scalings with trap frequency (u;" 1 ) [23] and mass (m" 1 ), these results are 



consistent with the present results for a size-scaling of d~ 4 . Another measurement has been 
made with 40 Ca + 0. For that experiment u m /2iT « 4 MHz and d ~ 700 /im and the heating 
rate was 0.005/ms. Compared to the present experiments and the Hg experiment, this is 
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also consistent with a d 4 scaling, although it is certainly unlikely that all systems had the 
same patch field environment. 



IV. CONCLUSIONS AND OUTLOOK 

We have measured heating from the ground state of trapped ions in different traps. The 
magnitude of heating rates and the results of the size-scaling measurements are inconsistent 
with thermal electronic noise as the source of the heating. The results do not indicate any 
strong dependence on trap-electrode material or on the type of trap potential (pondermotive 
or static). The rf voltage applied to the electrodes may play a role in heating, in as much 
as it may have an influence on patch-potentials. 

Since we have not identified the mechanism for the observed heating, it is difficult to say 
what path should be taken to correct it. If fluctuating patch-potentials on the surface of 
the electrodes are the cause, then further cleaning may be appropriate. Additionally, better 
masking of the the trap electrodes from the Be source ovens may help. 

The results coupled with those of other experiments |IU|,|5| , strongly indicate that bigger 



traps have smaller heating rates. This is not a surprise, but the strength of the scaling 
may be. With little sacrifice in the trap secular frequency (which ultimately determines the 
fastest rate of coherent manipulation) a dramatic decrease in the heating rate vs. logic gate 
speed appears possible using larger traps. 

We acknowledge support from the U. S. National Security Agency, Office of Naval Re- 
search and Army Research Office. We thank Chris Langer, Pin Chen, and Mike Lombardi 
for critical readings of the manuscript. 



V. APPENDIX: THERMAL ELECTRIC FIELDS 

We are interested in the thermal electric field power spectral density SEi{r,u>) generated 
from a specified volume of conductor. The conductor can be decomposed into a web of 
resistors each carrying current spectral density Si i = AksT / Ri (where we assume ksT ^> 
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huj). The resistance along the ith direction of an infinitesimal volume element is Ri = 
dl/(adA), where a is the conductivity, dl is the length along i and dA is the cross-sectional 
area. A Fourier component of current Ii(ui) through the volume dV = dldA gives rise to an 
electric dipole P%{oS) = Ii(u)dl/cj, thus the equivalent spectral density of electric dipole of 
the infinitesimal resistor is isotropic: Sp.(u)) = 4ksTadV/u 2 . 

The electric field from an electric dipole P(r', uj) oscillating at frequency u and position 



r' is 



Ei(r,u)= £ P,(r\uj)G l3 (r,r',u). (12) 

In this expression, Gij(r,r' ,u) is a Green function matrix, representing the ith component 
of electric field at position r due to the jth component of a point dipole at r' which satisfies 
the appropriate boundary conditions of the geometry. The electric field spectral density at 
position r is an integral over the dipoles in the conductor volume: 

Ah T r 

S Et (r,u) = -y a{r') £ \G tJ (r, r', u)\ 2 dV. (13) 

j=x,y,z 

The Green function satisfies Gij(r,r',u) = Gij(r',r,u), so the above integral can be 
interpreted as the Ohmic power absorbed by the conductor from the electric fields generated 
by a point dipole at position r. By energy conservation, this must be equivalent to the time- 
averaged power dissipated by a point dipole at r, which is related to the imaginary part of 
the Green function matrix Gij(r, r, uj) [24j]. This simplifies Eq. (|13D, leaving the fluctuation- 



dissipation theorem 



S Ei (r,u) = '^ Yl SmGijfar,*). (14) 



Agarwal solved Maxwell's equations for Gij(r', r, uj) for the simple geometry of an infinite 
sheet of conductor filling the space z < with the conductor-vacuum interface in the z = 
plane [25|| . Although this idealized geometry is far from any real ion trap electrode structure, 
rough scalings of the thermal fields can be relevant to real ion trap geometries. From Ref. 
the Green function matrix for this problem is diagonal with axial (z) and radial (p) 



components 
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G xz {z 1 z,u) = Gf™(u)+i 



q /w e — w 
o w V w e + w 



dq 



(15) 



G pp (z,z,u) 



Qfr 



w 



+ 



it; — 



dq, 



(16) 



In the above expressions, s(u) = Eq + ia /e^uj is the dielectric function of the conductor 
(in the low frequency limit), k = u/c, and wavevectors Wq and w (generally complex) are 
defined by w\ = k 2 — q 2 and w 2 = k 2 e — q 2 with wq > and Sm w > 0. The free space 
Green's function G^ ree (uj) has imaginary part Qtb G^ ree {uo) = k 3 /Qtteq and gives rise to the 
isotropic free space blackbody electric field fluctuations when substituted into Eq. (|T4|). 

The above integrals are significantly simplified in the "quasi-static" limit, where C 1 



and the conductivity is sufficiently high so that kS *C 1, where 5 = ^2c 2 Eo/utj is the skin- 
depth of the conductor. Despite these conditions, no restriction is placed on the value of 
z/5. We break the above integrals into two pieces. The first piece / fc has q < k with wq real. 
In the quasi-static limit, this piece can be shown to cancel the free space contribution to 
the transverse Green function Sm G pp (z, z, u>) while doubling the free space contribution to 
the axial Green function 3m G 22 (z, Physically, the presence of the conductor negates 

the transverse free space blackbody field while it doubles the axial blackbody field due to a 
near-perfect reflection. The second piece of the integrals f£° has q > k with wq imaginary. 
These pieces of the integral can be solved to lowest order in kz and k5. Combining terms 
and substituting the results into Eq. (|i~4D, the thermal electric field spectral density is 



2k B Tu 2 k B T 
3nEoC 3 Atcctz 3 \ 



1 1 Z A 

2 + Vi + ^ 



(17) 



Se d (z,uj) 



k n T 



8naz 3 \ 2 + V 4 + 5 4 ' 



(18) 



These expressions show that the thermal electric field noise scales as 1/z 3 for z -C 5 fll4 



but scales as 1/z 2 for z ^> 5 [§,01- At large distances z > ^JS/k (with kz <C 1), the axial 
field noise settles toward twice the free space blackbody value while the radial field vanishes. 
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This result is also reported in Ref. p6[ . The behavior is shown in Fig. [7], where Eqs. (0J 
and (|18j) have been substituted into Eq. |], giving the expected thermal heating rate for a 
9 Be ion trapped with molybdenum electrodes at room temperature. Note that the predicted 
heating rate at trap sizes typical in our experiments is significantly slower than the 0.1-1 
quanta/s rate predicted in pyn]- This difference comes from the choice of the value of the 
resistance in Eq. ||, which was chosen in |p|,[iT| as an absolute upper limit. 

When interpreting these results, only the rough scaling should be considered. Realistic 
ion trap electrode geometries are more complicated than a single infinite conducting plane, 
involving a more closed electrode structure. This generally requires a full numerical solution 
to the relevant boundary value problem. Moreover, we are usually interested in the electric 
field fluctuations at the center of the trap, where these fluctuations will be substantially 
different from those above an infinite plate. 
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FIG. 1. An example of heating rate data. The main graph shows n (= (n)) vs t, the delay 
between cooling and probing. The insets show Raman spectra from which n is extracted, according 



to Eq. 11. For the insets, P| is the probability that the ion remains in the | J.) state after application 
of a Raman probe of fixed duration with difference frequency Acj; rsb: red motional sideband, bsb: 
blue motional sideband. The sidebands shown are the 1st sidebands. The data are for trap 5 from 
Table | at 5 MHz secular frequency and n = 12 ± 2/ms. 
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FIG. 2. Schematic diagram of the electrodes of trap 3 (from Table |). The distance between 
traps 3a and 3b is 1.7 mm. Not to scale. 



19 



0.2 mm 





FIG. 3. Schematic diagram of the electrodes of the linear traps (traps 4, 5 and 6 from Table |). 
The traps are formed by evaporating gold (approximately 0.75 /mi thick) on an alumina substrate. 
The outer segmented electrodes are the endcaps, while the long unbroken electrodes carry rf. The 
axial direction (labeled z) is parallel to the rf electrode. The two separate trap wafers are spaced 
by 200 /mi for traps 4 and 5 and 280 /mi for trap 6 (spacers not shown). Schematic diagram not 
to scale. 
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FIG. 5. Heating rates vs. trap secular frequency in a) the elliptical ring trap 2 and in the micro 
linear-trap 4 and in b) linear traps 5 and 6. The only intended difference between traps 5 and 6 is 
the size. In all four data sets, the secular frequency was varied by changing a static potential only. 
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FIG. 6. Data from trap 3, showing heating rates vs. trap secular frequency, (a) Data set # 
1. The two points on the small trap at u> m /2ir = 5.3 MHz were taken with Raman cooling to 
n(t = 0) ~ and with Doppler cooling only to n(t = 0) ~ 2. Note that they give comparable 
results, as they should. The dashed lines show a w" 1 scaling, (b) Data set # 2. The small trap 
data were taken with an rf voltage of ~400 V and the big trap with ~600 V. The secular frequency 
was changed by tuning the DC bias between fork and ring electrodes. 
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FIG. 7. Expected heating rate vs. distance from electrode for thermal electronic noise, from Eqs. 
(|l7|), ( |i"8| ) and The numerical parameters are those for 9 Be + with molybdenum electrodes at 
10 MHz secular frequency and room temperature. 
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linear 


Au 
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1.1 
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TABLE I. Summary of traps. The size column is approximately the distance between the ion 
and the nearest electrode surface. £It is the trap rf-drive frequency. The heating rate (n) is for 
a trap secular frequency of 10 MHz, which in the case of traps 3a and 3b had to be extrapolated 
from data at lower trap secular frequencies. The two numbers quoted for 3a and 3b are for two 
different versions of the trap. See the text for a further discussion of the parameters. 
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